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Abstract 

Analytical solutions are obtained to the nonlinear equations describing 
neutrino oscillations when explicit neutrino-antineutrino asymmetries 
are present. Such a system occurs in the early Universe if neutrinos 
have a non-zero chemical potential. Solutions to the equations lead 
to a new type of coherent behavior governed by two modes. These 
bimodal solutions provide new insights into dense neutrino gases and 
into neutrino oscillations in the early Universe, thereby allowing one 
to surmise the flavor behavior of neutrinos with a non-zero chemical 
potential. 
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I. Introduction 



In this work we undertake a theoretical study of flavor oscillations in dense neu- 
trino gases. Such gases appear in physically interesting systems. An example occurs 
during the early Universe when neutrinos are self-interacting and fill space densely 

3. Likewise, during the final collapse of a supernova, neutrinos are emitted co- 
piously [0, 3. If neutrinos have masses and mix then oscillations among neutrino 
flavors can effect the physics of these systems. 

The treatment of neutrino oscillations in a dense gas is not so straightforward. 
The flavor behavior of a particular neutrino depends on the flavor content of back- 
ground neutrinos. However, the background neutrinos also oscillate. To know the 
flavor content of background neutrinos, it is necessary to know the oscillations of 
all individual neutrinos. Nonetheless, a Hartree-Fock-like self-consistent formalism 
has appeared which can handle the system @]. The behavior of the gas depends on 
certain statistical properties such as the energy distribution, the nature of density 
perturbations and the initial production of neutrinos. For the case in which only 
an energy distribution is involved, the formalism has been further developed for the 
neutrino-antineutrino gas in refs. 0, ||, 0] . The physical effects of neutrino oscillations 
on supernovae is treated in ref. p[. 

Since the system is self-interacting and nonlinear, one usually has to resort to 
numerical methods to determined the physics and flavor behavior. This approach 



was used in refs. P, 0; B 10) analyze neutrino oscillations in the early Universe 



for the case in which neutrinos have a chemical potential /ij, which is zero. In a gas 
for which /x^ = 0, the total number of antineutrinos is equal to the total number of 
neutrinos. It is unknown whether fi,y = in the early Universe. Indeed, the excess 
of electrons over positrons implies that the chemical potential for charged leptons is 
nonzero. Because this excess is tiny, being related to the baryon asymmetry of the 
Universe, the chemical potential for electrons is quite small. Nevertheless, because 
neutrinos are so dense at one second after the Big Bang, when the temperature of 
the Universe is about an MeV, even a tiny /i^ can have an effect. One purpose of 
the current work is to surmise the flavor behavior of neutrino oscillations in the early 
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Universe when /j.^, 0. Various aspects of neutrino oscillations in the early Universe 
are also addressed in refs. [12-20]. 

Thus one is led to consider the system when the total number of antineutrinos is 
not necessarily equal to the total number of neutrinos. We assume the gas is homo- 
geneous and isotropic. Under this assumption, the averaged spatial neutrino currents 
are zero and one needs only to consider neutrino densities. We also assume that the 
energy distribution is the most important statistical property of the gas. For simplic- 
ity, two flavors, electron and muon, are treated. In many regions of parameter space, 
the three-flavor case is accurately approximated by the two-flavor case. Finally, we 
consider situations in which hard scattering processes are much smaller than forward- 
scattering phase effects so that hard scattering can be and is ignored. This occurs 
for times greater than one second after the Big Bang and outside the neutrino sphere 
of a supernova. In these situations, GpE'^ -C 1, where Gp — 1-17 x 10~^^ MeV~^ is 
the Fermi coupling constant and E is the energy of a typical neutrino. The analytic 
results as well as the graphs displayed in the figures hold only for the case in which 
non-forward scattering can be neglected. 

Several parameters describe the system. Two important ones are A, which is the 
mass squared difference between the two mass-eigenstate neutrinos, and 9, which is 
the vacuum mixing angle: 

A ^ {ml -ml) , (1.1) 

and 

J^l = ^eL cos 9 — U^L sin 9 , 1^2 = I'eL siu 9 + U^L COS 9 . (1.2) 

Here, ui and 1/2 are the mass-eigenstate neutrinos, which have masses mi and m2, and 
UeL and are the left-handed electron and muon neutrinos, which are the states 
entering in the currents of the weak interactions. 

Let Ni, and Np be respectively the number of neutrinos and the number of an- 
tineutrinos in the gas. We use labels j and k to enumerate respectively the neutrinos 
and the antineutrinos. Thus the index j runs from 1 to Ni, and the index k runs from 
1 to Np. The energy of the jth neutrino and the kth antineutrino is denoted by 
and E . The system is put in a box of volume V. The neutrino and antineutrino 
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densities and rip are given by riy = Ny/V and np = Np/V. 

The wave function for jth neutrino is governed by a two-component vector in 
flavor space, i.e., 

wliere (respectively, t'^*^'^) is the probability that the jth neutrino is an electron 
neutrino (respectively, muon neutrino). Likewise, for the kth antineutrino, the flavor 
wave function is 

Since, under the assumption of no hard scattering, neutrinos and antineutrinos are 
not created or destroyed but merely change flavor, conservation of probability implies 
that 

^>^+« = l, i^f^^e' + = 1 ■ (1-5) 

The equations governing neutrino oscillations are easier to understand intuitively 

■k 



and visually by using vectors and w associated with Eqs. (|1.3| ) and (|1.4| ) given by 

m 



^ (z/fz/,^ - 2Re(z^^z/J), 2Im(z^fz7;)) , (1.6) 

where Re and Im denote the real and imaginary parts of a complex number. 

The equations governing the flavor behavior of a gas of self-interacting neutrinos 
and antineutrinos are ^ H, ^ 

^^^.Bi, ^^.-x5^ . (1.7) 

where Bl and are given by 

^^ = ^-K., Bt = ^,+V:. , (1.8) 



with 

A = A(cos2^, -sin2^,0) . (1.9) 

Here, 



3 



implements the interactions of neutrinos and antineutrinos among themselves. An 
asterisk on a vector indicates a change in sign of the third component so that w* = 
{w,,W2,-W3). In Eq. ([TTq ), 

(t;) = E^^ H = E^' • (1-11) 

j k 

The neutrino-neutrino term V^^, makes the system of equations in Eq. ( p..?] ) nonlinear. 

The equations in Eq. (|1.7|) resemble the motion of charged particles in "magnetic 
fields" Bi and 5^. 

The conservation of individual neutrinos in Eq. (|L^) corresponds to 

v^-v^{t) = l, w''-w''{t) = l . (1.12) 

The total number of neutrinos Ni, and antineutrinos Np is thus N^, = J2j and 
N, = Ek\w''\. 

If, at time t = 0, the jth neutrino and kth antineutrino begin as an electron 
neutrino and an electron antineutrino then 

t;^(0) = (1,0,0) , tZ;'=(0) = (1,0,0) . (1.13) 

If they begin as a muon neutrino and muon antineutrino then tP (0) = (—1, 0, 0) and 
w^{0) = (—1, 0, 0). If there are both electron and muon neutrinos initially, then only 
the excess is relevant for neutrino oscillations. Thus if the there are more electron 



neutrinos than muon neutrinos at t = then Eq. ( 1.13 ) is general. In the current 



work, the initial conditions in Eq. ( |1.13|) are referred to as initial flavor eigenstates. 

II. Reduction of the Problem 

The ratio of the strength of the vacuum term to the maximum possible neutrino- 
interaction term for the jth neutrino and kth antineutrino is given by the parameters 
and R^, where 

f.j = ^ j.k ^ _^ (2 1) 

V2GfE^ {n, + rip) ' ^GpE^ {n, + np) ' 

When K? <^ 1 and ^ 1, the vacuum term is much smaller than the nonlinear 
neutrino term. This corresponds to the case when the neutrino gas is dense. Under 
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these circumstances, ref. found that ahgnment holds. Ahgnment has been observed 
numerically 0] and is understood theoretically [^. When alignment holds, 

v^{t) ^ j| ^ m , w'{t) ^ M ^ f^^t) (2.2) 

are good approximations. Throughout the rest of this work, the dense neutrino case 
is treated, so that Eq. (|2.2|) hold quite accurately. 

When Eq. (|2.2| ) is an exact equality, the system of 3Ni, + 3Np differential equations 
in Eq. ( p..7|) reduce to six equations. One sums over j and k in Eq. (L7) to obtain 



where 

= V2Gp {n,f, - n,f:) , (2.4) 

and 1/Eq and 1/Eq are the average inverse neutrino and antineutrino energies given 
by 

Eo N,^E^ ' Eo N.^e" 
Initial flavor eigenstates correspond to 

r,(0) = (1,0,0) , f^(0) = (1,0,0) . (2.6) 

III. The Solution 

It is convenient to change variables. First, we use r in lieu of t, where 

T = -—j—t 3.1 

2Eo^Ko ^ ' 

with 

Ko = • 3.2 

Second, we use a mass eigenstate basis. This involves a rotation of vectors by an 
angle of 1Q. We denote the new vectors for neutrinos by y and for antineutrinos by 
z. They are given by 

?/i(r) = r^i(t) cos 26* — rt,2(t) sin 26* , 
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1/2 (r) =r^i{t) sin 26 + r^2 (t) cos 26 , 

y;{r) = r,3{t) , (3.3) 

and 



Tl 

zi{t) = — {ryji{t) cos 261 - r^2(i) sin26') 
n,- 



Tl 

Z2{t) = — (r^i(t) sin 261 + r^2W cos26') 
rip 



71 

z,{t) ^ -^r^.it) . (3.4) 
rip 



For antineutrinos, a factor of ny/rip is included for convenience, and the sign of the 
third component is flipped, which corresponds to using the "*" vector. Equation (|2.3| ) 
becomes 

^ = y X ^761 - i {y- z)\ , ^ = -zx Uei + ^ (y - ^ , (3.5) 



where 



7 = \f^o , 7 = ^7 • (3.6) 



Initial flavor eigenstates correspond to 



y(0) = (cos 26, sin 26, 0) , z(0) = — (cos 26, sin 2^, 0) . (3.7) 

rip 

Equation ( p.5|) is solved for arbitrary initial conditions in Appendix A. The solu- 
tion for yi and Wi is in terms of Jacobi elliptic functions. Thus the motion of the first 
components of y and w is period in r with a period which we denote by Tf. See Eqs. 
( |A.18| ) and (|A.2Q|) . The solutions for the second and third components are 



y2{T) = \/l-yi{r) cos ip{T) , y^ir) = y^l -?/f(r) sin^(r) , (3.8) 
Z2{t) = \/(~) -ih-yi{r)f cos^(r) , 



n ^ ^ 



n,- 



Z3{t) = J -) -{h-y,{T)y sin ^{t) , (3.9) 



where ki = yi (0) —zi (0) and where and ip are complicated functions of yi. Because 
yi is periodic, cos ip^r) and sin(y9(r) are periodic with a period 7^. Likewise, cos ipi^r) 
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and simplr) are periodic with a period See Eq. (|A.25|) . The periods in t are 
related to those in r by the factor 2Eq^/k^I /S. according to Eq. 

Hence, the motion of the second and third components is bimodaL They involve 
two periodic functions each with a different period. In terms of the original variables 
f\, and fv,, all components are bimodal because r^i and r^i involve ?/2 and W2 when 
Eqs. ( |3.^j| ) and ( |3.4| ) are inverted. In general, the periods are not compatible and the 
curves for r^i, r^,2, ^^3, r^i, r^2 and r^s appear irregular. When the gas has only 
neutrinos and no antineutrinos, the slow mode goes away and the behavior is purely 
periodic. This was the case observed in ref. [Q. Likewise, in the CP-symmetric case, 
in which = Np and Eq = Eq, the periods Tg and Tg go to infinity and the motion 
also becomes periodic with period 7}. This is the case treated in ref. [pT[ . 



IV. Explicit Examples 

The analytic solution in Appendix A is quite complicated. To obtain insight, it 
is useful to consider specific cases. Below, the values for sin 26, A, 

•} ^U"} etc., are 

not chosen to correspond to those in physical systems, such as in the early Universe, 
but to best illustrate different behaviors. 

It is convenient to use r in Eq. (|3.1| ) since then the system depends on only four 
parameters, namely, sin 26', 7, 7 and n^/rip. Although solutions also depend on the 
initial values of y (0) and -2(0), in this section, we avoid this dependence by using the 
initial flavor eigenstates in Eq. ( ^TTf ). We select sin 26' = 0.8. The large value for the 
vacuum mixing angle amplifies effects. 

Figure 1 shows the three components of y (t) for the case n^/rip = 0.9, 7 = 0.1 and 
7 = 0.12. Since the solution for yi (r) is always in terms of a Jacobi elliptic function, 
the graph of yi (r) is periodic. This is evident in Figure la. Since two periods are 
involved in y2 (r) and y^ (r) , the corresponding curves are bimodal and irregular, as 
Figures lb and Ic show. The plots for zi (r) , Z2 (r) and (r) are qualitatively similar 
to yi (r), y2 (r) and y^ (r). In other words, Zi (r) is periodic, and Z2 (r) and z^ (r) are 
bimodal and irregular. To save space, we do not display them. Additional insight is 
gained by examining the orbit for neutrinos, i.e., the trajectory {yi (r) , y2 (r) , 1/3 (r)) 
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in three-space. Figure 2a displays the neutrino orbit. The motion is as though one 
were moving back and forth over a cup-shaped surface in which the cup is also slowly 
rotating. These two motions crudely represent the two modes: The back-and-forth 
mode has a period of 7}. The other mode corresponds to the rotation of the cup, is 
governed by the period 7^, and is controlled by ip (r). The antineutrino orbit, given 
by -2 (t), is similar to yij) except that "little loops" are made near the "edge" of the 
cup. See Figure 2b. In the flavor eigenstate basis, all components are bimodal and 
irregular. Figure 3 displays r^i (r) = cos26'yi (r) + sm26y2 ij). Qualitatively r^n (r) 
is similar to r^i (r). The case in this paragraph is near the CP-symmetric limit for 
which n^jriy = 1.0 and 7 = 7. In that limit, the motion is back-and-forth over the 



cup but the cup does not rotate so that the motion is a single mode and periodic 
in all variables. 

When n^/riy is reduced, the amplitude for the back-and-forth motion becomes 
smaller and does not pass over the "top" of the cup. Figure 4a shows the orbit of 
y(r) for the case Uyjny = 0.75, 7 = 0.1 and 7 = 0.12. The antineutrino orbit is 
similar to the neutrino orbit and again has little loops. See Figure 4b. In terms of 
components, yi (r) and Zi (r) have graphs similar to the one in Figure la but the 
amplitude is smaller. The graphs for y2 (r), y^ (r), Z2 (r), (r) look like cosine and 
sine functions with "small bumps" . Figure 5, the plot of ?/2 versus r, exemplifies this. 

When riy/up = 0, one obtains the pure neutrino gas with no antineutrinos. The 
orbit for the neutrinos in this limiting case, which was analytically obtained in ref. 
2l|], is a circle. Figure 6a displays the neutrino orbit for the case n^/np = 0.4, 



7 = 0.5 and 7 = 0.6. As expected, the orbit is approximately a circle. However, 
because n^/np 7^ 0, there are "wiggles". The orbit looks somewhat like a "lasso". 
The antineutrino orbit is "flowerlike" . See Figure 6b. 

It turns out that, in the above three examples, % = Ts. For 7 and 7 not too large, 
this equality is general and holds for other values of the parameters sin 2^, nu/np, 
etc.. Hence the angular motions, (t) and <^(t), do not get out of phase. As a 
consequence, the initial CP-asymmetry is not amplified in time. For the discussion 
in the Section VI, we refer to this phenomenon as ip-phase locking. 

However, when 7 and 7 become sufficiently large, the equality 7^ = 7^ is suddenly 
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violated. For initial flavor conditions, this happens at the point in parameter space 
when 

Tt 

^(l-cos2e)=y+ , (4.1) 
rtr, 



where ?/+ is given in Eq. (|A.17 ). This is the point at which the denominator factor 



{ny/riyf — [yi — kif' in ^ vanishes for the maximum value obtained by yi. See Eq. 
( |A.23| ). For 7 or 7 above the critical point, Tg ^ Tg. In this case, even a small initial 
CP-asymmetry can eventually grow into a large CP-asymmetry at certain times. 

As an example of this phenomenon, consider the case with Uv/np = 0.9, 7 = 0.5 
and 7 = 0.5. One flnds Tf ^ 2.2832, Z ^ 4.8383, % ^ 4.3235. Figures 7a and 7b 
compare y2 (r) and Z2 (t). By the time r is 10, y2 and Z2 are about one-sixth of a 
period out of phase so that a sizeable CP-asymmetry is built up, even though the 
initial CP-asymmetry, as measured by the difference of n^/up from one, is 10 %. The 
orbits for this case are fairly irregular. See Figures 8a and 8b. 

V. Relation to Behaviors Observed in the Early Universe 

In a certain region of parameter space for A < 0, ref. |^ observed irregular flavor 
behavior. It is interesting to ask whether this could be bimodal self-maintained 
coherence. The equations governing neutrino oscillations in the early Universe 
1^ differ from those in Eq. ( |1.7| ) in two repects: First, there are CP-asymmetric 
effects. This means that the medium-induced linear term, the analog of A in Eq. 
( p..8| ) is slightly different for neutrinos and antineutrinos. Second, Bi and are 
time dependent due to the expansion of the Universe. However, in much of parameter 
space, interaction rates are faster than the expansion rate of the Universe, so that for 
short times B-l and B^ can be considered constant. Over long times, however, effects 
due to a non-constant Hubble constant can be important. 

The CP-asymmetric effects, which are due to the excess of electrons over positrons, 
are relatively small. However, they can eventually lead to a situation in which an- 
tineutrino vectors are not exactly aligned with neutrino vectors. Since the neutrino 
gas is extremely dense, small CP-asymmetries in the neutrino sector may lead to a 
sizeable interaction in V^,^. Hence, the electron-positron asymmetry might be impor- 
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tant indirectly via the nonlinear neutrino-neutrino term. 

Given this idea, it is reasonable to assume that the simplified model in Sect. II 
can represent the coherent behaviors for the early Universe, for the case where the 
energy distributions of neutrinos and antineutrinos are the same. The main effect 
of the positron-electron asymmetry is to misalign y and z. Hence, it is of interest 
to examine solutions to Eq. ( |2.3| ) with Eq = Eq, Uu = Up and |r„| = |fv,|, but with 
r.(0) ^r-;(0). 

We perform the same change of variables as in Eqs. (|3.3|) and (|3.4|) . The CP- 
asymmetric effects are then represented as a change in the initial conditions for z 
via 

z(0) = — (cos (20 -2a), sin (20 -2a) COS 2/5, sin (20 -2a) sin 2/5) , (5.1) 



while keeping y(0) as in Eq. ( p.7|) , so that y(0) = (cos 20, sin 20, 0). Here, a and [3 
are angles that control the CP-asymmetry in the neutrino sector in the 1-2 plane and 
in the 3 direction. 

Figure 9 displays the orbits for the neutrinos and antineutrinos for sin 20 = 0.8, 
sin 2a = 0.01, sin 2/5 = 0.0, 7 = 7 = 0.1 and nyjuy = 1.0. The motion is mostly 



back-and-forth as in the CP-symmetric case of reference |21| but the second mode 
causes the back-and-forth cycle to shift somewhat and not retrace over itself. When 
both a and (3 are non-zero, fairly irregular behavior can arise. See Figure 10 which 
displays the orbits for the neutrinos and antineutrinos for sin 20 = 0.8, sin 2a = 0.1, 
sin 2/5 = 0.1, and 7 = 7 = 0.1. These examples show that bimodal self-maintained 
coherence and irregular fiavor dependence arise even when the only CP-asymmetry 
is in the initial conditions for neutrinos and antineutrinos. 

We have examined some of the orbits in the simulations of refs. |Tl| to see 
whether bimodal coherence occurs in the early Universe. Due to the variety of mani- 
festations of bimodal solutions, as typified in Figures 1-10, it is not always so easy to 
be sure that the irregular behavior seen corresponds to bimodal coherence. However, 
many results from oscillations in the early Universe are similar to the examples consid- 
ered here. One clear case from ref. |ll| occurs for the simulations with A = — l.Ox 10^^ 
eV^, sin^ 20 = 0.49 at approximately 0.13 seconds after the Big Bang. Figure 11 dis- 
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plays the orbit for the neutrinos in the flavor basis. The orbit for antineutrinos is 
virtually the same because of the CP-suppression mechanism discovered in refs. ||^, 0. 
Two modes are clearly seen. The orbit in Figure 11 is qualitatively similar to the 
orbit in Figure 2. Another example from ref. [Tl|] occurs for A = —1.0 x 10^^ eV^, 



sin^2^ = 1.0 X 10^*^ at approximately 0.94 seconds. Although we do not display the 
orbit, it is similar to Figure 9a. 

Bimodal self-maintained coherence might be the general behavior of many dense 
neutrino gas systems. All three behaviors observed in the simulations of refs. ^ |ll| 
are qualitatively obtained as limits. Irregular behavior corresponds to bimodal coher- 
ence when both modes have sizeable amplitudes. Residual self-maintained coherence 
occurs when only one mode has a sizeable amplitude. Finally, smooth behavior arises 
when vectors point in the directions of the the linear terms in the effective magnetic 
fields, so that both modes have small amplitudes. Indeed, for a dense neutrino gas 
in the early Universe, tiny rapid fluctuations were observed in the simulations of refs. 

^ ^, |10], |lT| in the smooth region of the phase diagrams. One possible explanation 
for them is that they correspond to small-amplitude bimodal coherence. Of course, 
when the amplitudes for oscillations are tiny, other effects may create them. Even if 
all the behaviors in parameter space for A < correspond to bimodal coherence, it is 
still useful to characterize them by irregular, residual self-maintained coherence, and 
smooth as in Figure 1 of ref. |p. 

VI. Neutrino Behavior in the Early Universe for ^ 



The studies of neutrino oscillations in the early Universe in refs. [^, |^, |^, |T0|, p 



were performed under the assumption that = Up and that the energy distributions 
for neutrinos and for antineutrinos are the same. Given that there is a CP-asymmetry 
for charged leptons, namely the excess of electrons over positrons, it is possible that 
the numbers of neutrinos and antineutrinos are not equal, in which case a chemical 
potential potential /i,^ for neutrinos needs to be introduced. 

One purpose of the current work is to gain insight into the asymmetric neutrino- 
antineutrino gas. With this insight, it is possible to speculate on the phase diagram 
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and on the flavor behavior. To rigorously determine these, simulations must be per- 
formed, but such computer studies are numerically intensive because one must explore 
a phase space of three parameters: 9, A and ^y. A thorough study cannot be under- 
taken with the power of present computers. However, an educated guess as to the 
results can be made. 

Due to the explicit asymmetry between neutrinos and antineutrinos when ^ 0, 
residual self-maintained coherence, which corresponds to a situation in which only 
a single mode is excited, is unlikely to appear, except for perhaps brief periods for 
exceptional values of parameters. Instead, irregular flavor, as governed by bimodal 
coherence, should replace it. Thus the diagram for A < probably has only two 
phases: smooth and irregular. When /x,^ 7^ 0, irregular behavior should arise in the 
regions of the phase diagrams of refs. |rT| where irregular or residual self-maintained 
coherence occurred. 

For A > 0, it is likely that both neutrino and antineutrino vectors follow their 
respectively effective magnetic fields, thereby maintaining themselves in approximate 
nonlinear mass eigenstates. Hence, smooth behavior is generally expected for A > 0. 
However, if Eq is significantly different from Eq then transitory weak irregular be- 
havior should take place in the "cross-over region". The cross-over region is the 
period during which neutrinos evolve from approximate flavor eigenstates to approxi- 
mate vacuum mass eigenstates 0. Small amplitude irregular behavior should appear 
there because neutrinos and antineutrinos make the transition to mass eigenstates 
at different times. This cause neutrino and antineutrinos vectors to point in differ- 
ent directions. The effect should be most pronounced for small A, particularly when 
A < 10^^ eV^. After the transition region, the behavior should become smooth again. 
A similar transitory effect might also occur in the smooth phase for A < 0. 

One motivation for considering a non-zero neutrino chemical potential is that it 
introduces additional CP-asymmetry into the system. Reference showed that a 
large neutrino-antineutrino asymmetry can have an effect on nucleosynthesis. How- 
ever, for the values of the parameters of the early Universe, <y9-phase locking is ex- 
pected. This implied that the CP-asymmetry due to ^ Q should not be amplified 
in time. Hence, neutrino oscillations for /Xj^ 7^ are likely to affect nucleosynthesis 
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at the same level as in the case of no neutrino masses and mixing |22, |2^. Con- 
straints on the chemical potential of neutrinos from Big Bang nucleosynthesis have 
been obtained in ref. p!6| . 



VII. Summary 

In this work, we analyzed neutrino oscillations in a dense gas in which there is an 
explicit neutrino-antineutrino asymmetry. The energy distribution for antineutrinos 
was not necessarily the same as the energy distribution for neutrinos, nor did the 
numbers of antineutrinos and neutrinos have to be the same. 

In the dense gas limit, alignment is expected to hold. Using this property, the non- 
linear multi-equation system was simplified to six equations in Section II. Remarkably, 
the six nonlinear equations were solvable. The solution was still quite complicated. 
After performing some transformations in Section III, most solution details were rel- 
egated to Appendix A. The generic behavior was bimodal self-maintained coherence. 
This is a collective soliton-like solution in which neutrino behavior is cooperative and 
in which two modes of different frequencies are present. 

Bimodal self-maintained coherence exhibits a wide range of behaviors depending 
on parameters and initial conditions. Particular illustrative cases were presented in 
Section IV. 

When the amplitude for one of the modes vanishes, ordinary self-maintained co- 
herence is obtained. Such behavior was first observed in test simulations in ref. [Q. 
It was subsequently observed in numerical simulations of the early Universe in ref. [H 



and understood analytically in ref. ||21[. Another neutrino flavor behavior, which has 
been seen in the early Universe, is smooth 00]. It corresponds to a bimodal solution 
in which the amplitudes of both modes are small. Finally, irregular behavior, which 
can occur for an inverted mass hierarchy [^], is mimicked by bimodal coherence when 
both amplitudes are present. These results suggest that bimodal coherence may be 
the general flavor behavior for many dense neutrino gases. As discussed in Section 
V, our results provide new insight into the simulations of neutrino oscillations in the 
early Universe performed in refs. [§|, 0; ^ 0' 0- 
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With this additional understanding, it was possible to surmise the neutrino flavor 
dependence in the early Universe when neutrinos have a chemical potential, so that 
there is an explicit neutrino-antineutrino asymmetry. The expected behavior for 
A > is smooth with a possible weak bimodal coherent phase for a short time 
period. The expected behaviors for the inverted mass hierarchy case, A < 0, are 
smooth and bimodal coherence. Bimodal coherence probably arises where irregular 
or residual self- maintained coherence arose for the fi^ = case of refs. [0, n|. Because 



of the phenomenon of y^-phase locking, effects on nucleosynthesis are not expected to 
be amplified by neutrino oscillations, as one might have a priori thought. The effects 
are probably limited to those directly due to a non-zero neutrino chemical potential 

In short, our work has shown that dense neutrino gases can exhibit interesting and 
non-intuitive phenomenon. The gases often have collective, cooperative and coherent 
flavor behavior. 
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Appendix A 

It turns out that Eq. ( ^.5| ) can be solved for arbitrary initial conditions. Four 
conserved quantities play a role. The first two 

l/l' + 2/2' + l/3' = l (A.l) 



2 



z,' + z,' + z,'= ^ (A.2) 



and 



correspond to conservation of neutrinos and of antineutrinos. Equation ( p.5|) also 
implies that 

yi- zi = ki , (A. 3) 
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and that 

y ~ 27 (7V 7) " "^'^^ + ~ '''^') " ' ^^-^^ 
where ki and ki2 are time-independent constants. They are determined by initial 
conditions, i.e., 

fci = yi(0)- 2:1(0) , 
= VM - 2^ ^) ((^2(0) - Z2(0))' + (2/3(0) - ^3(0))') . (A.5) 
For initial flavor eigenstates, 

/ ft 
ki cos 26 [1 



sir? 29 /, n,, \ ^ 



fci2^cos20 1-— . (A.6) 

27(7 + 7)V nj ^ ' 

It turns out to be useful to solve for the quantity ?/2 ^2 +2/3-23 using Eqs. (|A.1| )- (|A.4|) . 
One finds 



2/2^2 + 2/32:3 = ^ ^ + (~) ) + 7 (7 + 7) {ki2 - 2/1) + kiyi - yi 



(A.7) 



Thus, Eqs. ( |A.1| )- (|X^ ) allow one to eliminate four combinations of variables in 
favor of 2/1: 

2 , 2 1 2 
2/2+2/3-^1-2/1 , 



zi-^yi-ki , 
^2^ + ^3^ - f-V - (2/1 - fci) 



nr. 



2/2-22 + 2/3-23 ^]^[l-ki^+ ^ + 7 (7 + 7) (fci2 - 2/l) + ^ll/l " l/i' • (A.8) 

Take the first order differential equation for 2/1 in Eq. ( |3.5| ), differentiate it with 
respect to r and use the equations for the second and third components in Eq. (|3.5|) . 
It is a miracle that precisely the four quantities on the left-hand side of Eq. ( |A.8| ) 
enter. Straightforward but lengthy algebra gives 

(fyi ( l\ , 

= I 1 + - I (2/22:2 + 2/32:3) + 
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= -^Ayl-2By^-C , (A.9) 

where, in the last equahty, the substitutions in Eq. ( |A.8D are used. The constants A, 
B and C are given by 

7 

A _ (7 + 7)' fcl' (7 + 7) (fcl + fcl2) 

2 ^272 7 
+ (2±jl{k,k,,_-\{l-k^+{'^']\\ . (A,10) 



7 V 2 \ Knti 

Thus the motion of yi corresponds to a classical particle moving in a one-dimensional 

cubic potential. This system can be solved by quadratures. Multiply Eq. ( |A.9D by 

^ and integrate to obtain 

\{^f)-viy.) . (A.n, 

where 

V (yi) = -\ (^(0)) +C{y,- yM) + ^ (yi^ - (yi(0))^) + i {y,^ - {yMY 

(A.12) 

The potential V (yi) has three zeroes and is negative between the two larger zeroes. 
The motion of yi is between these two larger zeroes where V {yi) < 0. Integration of 
Eq. (^[np leads to 

r=r'^'-jL= . (A.13) 
Jy^io) ^-2V{y) 

This equation gives r in terms of ?/i(t) where the latter appears as the upper limit of 
integration. Equation ([A.13| ) is inverted to obtain yi as a function of r. 

Since the potential is cubic, the solution is always related to Jacobi elliptic func- 
tions. In general, a cubic equation must be solved, so that the explicit solution is 
rather complicated. However, when 

^(0) = , (A.14) 
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one can avoid this difficulty. The initial condition in Eq. ( A.14|) is an important case 
because it includes the initial flavor eigenstates in Eq. (p.7|). When Eq. ( [A.14| ) holds, 
the potential becomes^ 

V {yi + vM) = yi {Ayl + By, + C) = (yi - y+) {y, + y^) . (A.15) 

where 

A = A , 
B = B + 3Ay,{0) , 

C = C + 2By,{0) + 3Aiy,{0)f . (A.16) 
The zeroes of the potential are yi{0) — ?/_, yi{0) and yi{0) + y+ where 



-B + VB^ - AAC B + - AAC , , 

= 2A ' y- = — 2A — ■ ^^-^'^ 

Thus yi oscillates between yi{0) and yi{0) +?/+. 
The integral in Eq. ( |A.13| ) is 



Vfir)=ym^'-^ \ . (A.18) 

where 

Vfi2 - AAC 

yo = J = y+ + y- , (A.i9) 

and where sn and dn are the sine-amplitude and delta-amplitude Jacobi elliptic func- 
tions. 

Hence yi undergoes periodic motion with a period 7} in r given by 



where 



dw 



K i-e) = , , (A,21 

^ Jo ^(1 - ti.2) (1 + qV) 



^ For the case of general initial conditions one must shift yi by a constant which renders the 
middle zero at 0. Then one proceeds as in the special case explicitly carried out here. 
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with = y+lu-- The period Tj in t is 

Tf = 2^^,Tf . (A.22) 

The solution for zi{t) is given by using Eqs. ([A .31 ) and ( [A.18| ). 
Finding the solution of the second and third components of y and z is not so 
straightforward. One proceeds by using the polar coordinate angles in the two-three 
plane, i.e., let ip = tan~^ (1/3/1/2) and (f = tan~^ (23/22). Differentiate ip with respect 
to r and use Eq. (|3.5|) . Do the same for ip. Miraculously, only the four quantities in 
Eq. ( |A.8|) enter. Hence, the substitutions in Eq. ( |A.8| ) can be used to express ^ and 
^ in terms of 2/1 (t). After straightforward but lengthy algebra, one finds 



dr 



d^__ fei+(7(7 + 7) (fci2-yi)-|(l + fci^-te)' 
dr ^ 7(l-?/i2) 



\ (1 - fci' + (^)') + 7 (7 + 7) (fci2 - yi)) ih - yi) + {^fyi 

(A.23) 



te) -(1/1-^1; 



Thus, 



^(r) = ^(0) + [ dr^ , ^(r') = ^(0) + dr^ , (A.24) 
where ^ and ^ are given in Eq. ( |A.23| ) and where 

The solutions for the second and third components are then given via Eqs. ( |3.8| ) and 

Since yi is periodic, the angular motion in the two-three plane is also periodic 
but with a period different from Tf. For neutrinos, let Tg denote the period in r. For 
antineutrinos, let % be the period. Then, 
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where | | indicates absolute value. The subscripts / and s on 7} and Tg stands for 
"fast" and "slow". In the extreme dense gas hmit 7} is much smaller than 7^, i.e., 7} 
is the period for the faster motion. The periods Tg and Tg in terms of t are 

Tg^2^^^%, Tg^2^^Tg . (A.26) 
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Figure Captions 

Figure 1. Components of y as a function of scaled time r for < r < 20 for the 
case sin 2^ = 0.8, ny/np = 0.9, 7 = 0.1 and 7 = 0.12 and with initial conditions as in 
Eq. ( p.7| ). (a) The component yi. (b) The component 1/2. (c) The component 1/3. 

Figure 2. The three-dimensional orbits for the example in Figure 1. (a) The 
neutrino orbit, (b) The antineutrino orbit. 

Figure 3. The flavor component function of r. Parameters are the same 

as in Figure 1. 

Figure 4. The orbits during < r < 10 for the case sin 2^ = 0.8, nj^/rip = 0.75, 
7 = 0.1 and 7 = 0.12 and with initial conditions as in Eq. (p.7|) . (a) The neutrino 
orbit, (b) The antineutrino orbit. 

Figure 5. A plot of 1/2 versus r for the example in Figure 4. 

Figure 6. The orbits during for < r < 20 for the case sin 2^ = 0.8, n,,/np = 0.4, 
7 = 0.5 and 7 = 0.6 and with initial conditions as in Eq. ( p.7|) . (a) The neutrino 
orbit, (b) The antineutrino orbit. 

Figure 7. A comparison of y2 and Z2 when yj-phase locking does not hold for the 
case sin26' = 0.8, n^/rip = 0.9 and 7 = 7 = 0.5. (a) The component 2/2. (b) The 
component Z2. 
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Figure 8. The orbits during < r < 10 for parameters as in Figure 7. (a) The 
neutrino orbit, (b) The antineutrino orbit. 

Figure 9. The orbits during < r < 10 for the case sin 20 = 0.8, n^jn^ = 1.0 
and 7 = 7 = 0.1 and using the initial conditions in Eq. ( ^.1|) with sin 2a = 0.01 and 
sin 2/5 = 0.0. (a) The neutrino orbit, (b) The antineutrino orbit. 

Figure 10. The orbits during < r < 10 for the case sin2^^ = 0.8, riu/rip = 1.0 and 
7 = 7 = 0.1 and using the initial conditions in Eq. ( p.l[ ) with sin 2a = sin 2/5 = 0.1. 
(a) The neutrino orbit, (b) The antineutrino orbit. 

Figure 11. A neutrino orbit with axes in the flavor basis for A = —1.0 x 10~^ eV^ 
and sin^ 20 = 1.0 x 10^^ during a small time interval at approximately 0.13 seconds 
after the Big Bang. 
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